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Abstract 

Background: Phenome-Wide Association Studies (PheWAS) can be used to investigate the association between 
single nucleotide polymorphisms (SNPs) and a wide spectrum of phenotypes. This is a complementary approach to 
Genome Wide Association studies (GWAS) that calculate the association between hundreds of thousands of SNPs 
and one or a limited range of phenotypes. The extensive exploration of the association between phenotypic 
structure and genotypic variation through PheWAS produces a set of complex and comprehensive results. Integral 
to fully inspecting, analysing, and interpreting PheWAS results is visualization of the data. 

Results: We have developed the software PheWAS-View for visually integrating PheWAS results, including 
information about the SNPs, relevant genes, phenotypes, and the interrelationships between phenotypes, that exist 
in PheWAS. As a result both the fine grain detail as well as the larger trends that exist within PheWAS results can 
be elucidated. 

Conclusions: PheWAS can be used to discover novel relationships between SNPs, phenotypes, and networks 
of interrelated phenotypes; identify pleiotropy; provide novel mechanistic insights; and foster hypothesis 
generation - and these results can be both explored and presented with PheWAS-View. PheWAS-View is freely 
available for non-commercial research institutions, for full details see http://ritchielab.psu.edu/ritchielab/software. 
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Background 

In Phenome-Wide Association Studies (PheWAS), the 
association between single nucleotide polymorphisms 
(SNPs) and an extensive range of phenotypic measure- 
ments are calculated in a high throughput, unbiased 
manner. The phenotypic data used in PheWAS can 
come from a variety of sources. One possible source is 
epidemiologic health surveys linked to genotypic data 
that include measurements of intermediate traits or bio- 
markers such as blood cell counts and blood pressure 
measurements, as well as information on case/control 
status for multiple clinical conditions and risk factors 
such as presence/absence of diabetes or hypertension. 
One such example is the Population Architecture Using 
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Genomics (PAGE) network, which is a National Human 
Genome Research Institute (NHGRI)-supported network 
of four study sites and a coordinating center accessing 
eight extensively characterized studies for PheWAS 
studies in diverse populations [1,2]. These survey-based 
PheWAS efforts are complimentary to on-going Phe- 
WAS efforts using electronic medical records linked to 
biorepositories such as those in the electronic Medical 
Records & Genomics (eMERGE) network [3,4]. 

The exploration of data in a PheWAS effort presents 
several challenges, including the need for data visuali- 
zation to assist with interpretation of the data. GWA 
studies of a single or limited number of traits lend them- 
selves to Manhattan plots where p-values for every 
test of association are plotted by chromosomal location 
(x-axis) and the level of significance is visualized easily 
(y-axis). Such a plot does not present the complex rela- 
tionships that exist between both genotypes and pheno- 
types in PheWAS. Therefore, to visualize the complex 
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Table 1 Listed are all the PheWAS-View plotting options, 
parameters, and flags (format: -flag name) for creating 
PheWAS-View plots using phewas_view.rb 

Usage: phewas_view.rb 
Standard PheWAS Plot 
-v, -version 
-e phewas file 



Table 1 Listed are all the PheWAS-View plotting options, 
parameters, and flags (format: -flag name) for creating 
PheWAS-View plots using phewas_view.rb (Continued) 



-o output name 
-t title 

-f image type 

-w, -lowres 

-a, -rotate 

-p p-value threshold 

-m, -maxp maximum p-value 

-R, -redline p-value 
-b 

-A, -samp-size 

-I ancestry map file 

-c phenotype class names 

-B, -showbest 

-x phenotype/SNP file 
-r ancestry 
-s SNP ID 

-L phenotype list file 
-N, -no-lines 

-C phenotype correlation file 

Sun Plot Settings 
-S 

-s SNP ID 



Show PheWAS-View version 

PheWAS-View formatted file 
for input 

Optional output name for 
the resultant plot 

Main title for the plot 
(enclose in quotes) 

Image format for output 
(png default). Other options 
depend on ImageMagick 
installation. 

Low resolution image 
(72 dpi) 

Rotate final image 
90 degrees 

p-value threshold, values 
less significant will be 
plotted in grey 

Maximum p-value to plot. 
Values less significant than 
the specified cut off are 
not plotted 

Draw a red line at the 
designated p-value 

Include direction of effect 
on plot 

Include sample size plot 

Optional ancestry map file 

Only results matching this 
phenotype class name are 
plotted 

Display detailed information 
for best score at each 
phenotype 

PheWAS expected SNP/ 
phenotype file 

List of race/ethnicities to 
include (AA, EA, MA) 

SNP ID to display from 
input file 

Optional phenotype list for 
inclusion 

No background lines drawn 
on plot 

Optional file with phenotype 
correlations 



Produce sun plot 

SNP ID to display in center 
of sun plot 



-P phenotype name 


Phenotype to display in 




center of sun plot 


-g gene name 


Gene to display in center 




of sun plot 


-G 


Include gene name along 




with SNP, when SNP is 




SclcLlcU IOI SUM piOl 


-E 


Include ancestry as 




description of result for 




sun plot 


-m, -maxp p-value 


Choose a p-value threshold, 




p-values less significant will 




not be plotted 


-p p-value 


For plotted results, any results 




more significant will be plotted 




in red 


-b 


To apply direction of effect for 




phenotypes in sun plot, - is 




negative direction, + is positive 




direction 



results of PheWAS, we have developed PheWAS-View, 
software that can be used to create visual summaries of 
the SNP, gene, phenotype, and association information 
resulting from these studies. Using specialized tools such 
PheWAS-View to investigate results on a larger summary 
level as well as the individual result level is key for inter- 
pretation, analysis, and sharing of PheWAS results. While 
this tool was developed specifically for PheWAS, it could 
be used in other high throughput bioinformatics data 
where thousands of association results are being explored. 

Implementation 

PheWAS-View was developed in Ruby, using the RMa- 
gick graphics library, for use at the command line. 
Through the use of various options, various output plots 
are possible. Table 1 shows various commands and op- 
tional settings for PheWAS-View. 

A single input file is required to produce a standard 
PheWAS-View plot (example Additional file 1). Required 
columns for the standard input file include a column of 
unique SNP identifiers such as an rs number, a column of a 
unique phenotype description/identification for the tests of 
association that were calculated for each SNP, and a column 
of p-values for each test of association. By adding additional 
columns, additional features are possible with PheWAS- 
View (example Additional file 2). Table 1 lists the various 
parameters/flag settings available for modifying PheWAS- 
View plots (format: -flag name) at the command line. 

Results and discussion 

One way to inspect initial PheWAS results is first through 
visualizing all association results across phenotypes using 
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Figure 1 Standard PheWAS-View Output Example Plot. A series of simulated PheWAS results plotted in PheWAS-View for a group of 
phenotypes. The y-axis presents -log 10 (p-value) of the tests of association for all SNPs for each phenotype, and the x-axis represents individ 
phenotypes. In this way all results for all SNPs are plotted for each phenotype, 



ual 



PheWAS-View. Figure 1 shows simulated PheWAS 
results plotted in standard PheWAS-View format for a 
series of phenotypes, using the example Additional file 1. 
The plot is similar in style to a Manhattan plot, where the 
y-axis represents the magnitude of the association results 
in -log 10 (p-value). However, unlike a Manhattan plot 
where the x-axis represents genetic location, the x-axis of 
a PheWAS-View plot represents each phenotype from 
the tests of association, plotted in the order phenotypes 
are listed in the PheWAS-View input file. This file can be 
sorted in a number of ways such as by ascending or des- 
cending p-value, and then plotted in PheWAS-View, 
(such as Additional file 3: Figure SI) to aid in exploration 
of the results. 

PheWAS-View options allow for modified views of 
results that may aid in further result investigation and 



interpretation. For instance, with an extensive number 
of phenotypes, the plot may become very large. One op- 
tion for managing plot size is using a filter based on a 
phenotype group or class. Phenotypes that are related 
can be given a unique group identifier. If this informa- 
tion is supplied in a column titled "phenotype_class", the 
PheWAS-View output can be filtered on any phenotype 
class of choice (using command line parameter -c pheno- 
type class). In this way, only results for phenotypes 
within that phenotype class are plotted. Figure 2 shows 
an example filtering on the phenotype class "Allergy" in 
PheWAS-View using Additional file 2, thus limiting all 
the plotted results to allergy related phenotypes. An 
alternate approach is to use the -L flag, and supply a list 
of phenotypes to filter the resultant PheWAS-View 
plot for, whereby the results for only those specific 
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Figure 2 Selecting Phenotypes for the Plot. PheWAS-View allows the user to filter the data on various criteria, including limiting the plot to a 
specific group or list of phenotypes. This is an example of filtering on the phenotype class "Allergy" in PheWAS-View. The data plotted only 
include those phenotypes that were notated as being in the allergy phenotype class. PheWAS-View has an option for supplying two pieces of 
phenotypic information, a short and a long phenotype description, in this case the short phenotypic description for all results is "Phen". 
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Figure 3 (See legend on next page.) 
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(See figure on previous page.) 

Figure 3 Options for Highlighting Results Based on p-value Thresholds. Figure 3A shows the result of using the parameter -p p-value, 
where results more significant than a specified p-value threshold are plotted in blue (in this case p = 0.01 ), and the other results are plotted in 
grey. Figure 3B shows the result of using the parameter -R p-volue to plot a red line at a p-value of interest (p = 0.01). Figure 3C shows the results 
of using the parameter -m p-value to plot only those values more significant than a chosen p-value threshold (p = 0.01). 



PheWAS 



-logltyp value) 



Phen: Alternaria flare lengtl 
Phen: Alternaria ■ flare widtl 
Phen: Alternaria • wheal lengtl 
Phen: Alternaria - wheal widtl 
Phen: Bermuda grass ■ flare lengtl 
Phen: Bermuda grass - flare wirJtfn 
Phen: Bermuda grass - wheal lengtl 
Phen; Bermuda grass wheal widtl 
Phen: Cat flare lengtl 
Phen: Cat - flare widtl 
Phen; Cat ■ wheal lengtl 
Phen: Cat - wheal widtl 
Phen: German cockroach - flare lengtl 
Phen: German cockroach ■ flare widtl 
Phen; German cockroach - wheal lengtl 
Phen: German cockroach - wheal widtl 
Phen: Mile - flare lengtl 
Phen; Mite - flare width 
Phen: Mile - wheal lengtl 
Phen: Mite - wheal width* 
Phen: Negative control ■ flare length* 
Phen: Negative control - flare width* 
Phen: Negative eonlrol - wheal length* 
Phen: Negative control ■ wheal width* 
Phen: Peanut - flare length* 
Phen: Peanut - flare width* 
Phen. Peanut - wheal length* 
Phen: Peanut ■ wheal width* 
Phen: Ragweed - flare length* 
Phen: Ragweed flare width* 
Phen; Ragweed ■ wheal length* 
Phen: Ragweed - wheal width* 
Phen: Russian thislle ■ ffare length* 
Phen: Russian I his He - flare width* 
Phen: Russian thislle wheal length* 
Phen: Russian thistle ■ wheal width* 
Phen: Rye grass flare length* 
Phen: Rye grass - flare width* 
Phen: Rye grass ■ wheal length* 
Phen: Rye grass - wheal width* 
Phen: While oak flare length* 
Phen: While oak - flare width* 
Phen: While oak - wheal length* 
Phen: While oak wheal width* 
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rs46O7103 ADAMTS9 - 
rs 10757276 CDKN2A/GDKN2B + 
rs 12596776 CETP 
rs43S580l TNRFSF11B <■ 
rs2363208 CDKN2A/CDKN2B + 
rs 10689353 ANG PTL3 , DOCK7 + 
rs3093058 CRP + 
rs2650000 HNF1A + 
rs6548238 TMEM18 - 
rs7679 FLIP + 

rs646?76 CELSR2 PSRClvSORT1/lpl3,3 - 

rs65ll720 LDLR 

rs150190S T1MD4 + 

rs3736228 LRP5 + 

rs 11 206510 PCSK9 - 

rs653689l LPL 

rs2237695 KCNOt + 

rs 13266634 SLC30A8 + 

rs 12654264 HMGCR + 

rs6S55911 GLUT9 ■ 
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rs7498665 SH2B1 * 

rs 17 145736 ML XI PL 
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rs560887 G6PC2,G6PC2 - 
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rsi2740374 CELSR2 - 

rs964184 APOA1.APOA1 ; APOA1/C3/A4/A5 

rs 1417938 CRP + 

rs460?103 ADAMTS9 - 

rs 141 7938 CRP + 

rs 17366743 ADIPOQ 

rs 1260326 GCKR.GCKR 

rs646776 CELSR2/PSRC1/SORT1/1 pl3.3 + 

rs3690162 ABCA1 



Figure 4 Vertical Format for PheWAS-View Plots. PheWAS-View allows results to be plotted in a vertical format. Figure 4 shows the same data 
used in Figure 2, plotted in a vertical format through using the parameter -a. Compared to previous figures in this manuscript, the phenotypes 
are now listed along the y-axis, and the x-axis represents -log 10 (p-value) of the tests of association. In addition, using the parameter -B will plot 
the SNP identifier as well as direction of genetic effect for the association (+ for positive direction of effect, - for negative) for the most significant 
association for each phenotype listed. Plotting PheWAS results this way facilitates the reading of phenotype descriptions, while still allowing the 
most significant SNP-phenotype results to be inspected visually. 
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Figure 5 Additional Information on Sample Size and Direction of Effect. For this plot the -b parameter was used to plot the direction of 
effect for all associations for each phenotype in a separate track. The sample size for each association was also plotted in a separate track using 
the -A parameter. 



phenotypes are plotted. Instead of (or in addition to) fil- 
tering results by phenotype group, data can also be fil- 
tered by SNP rsID using the parameter -s SNP ID. 

Multiple PheWAS -View options exist for highlighting 
results more significant than a specific p-value threshold, 
which may highlight results of interest. Using the param- 
eter -p p-value, results that are more significant than a 
specific threshold are plotted in blue, and the other results 
are plotted in grey (Figure 3A). Alternately, a red line can 
be applied at a p-value of interest (-R p-value) (Figure 3B), 
or through using (-m p-value) to plot only results more 
significant than a p-value threshold (Figure 3C). 

A useful alternate view is to plot the same information 
in a vertical format (-a), where the phenotypes are listed 



along the y-axis, and the x-axis are the -log 10 (p-value) 
from the tests of association (Figure 4). In this format, 
reading the phenotype identification is easier, while the 
most significant SNP-phenotype results can still be iden- 
tified visually. 

For plots in vertical format, using the parameter -B 
will plot the SNP identifier, gene symbol, as well as dir- 
ection of the genetic effect (positive (+) or negative (-)) 
for the most significant p-value for each phenotype 
(Figure 4). To plot effect size, these data must be provided 
in a column "ES", and gene symbol must be provided in a 
column "Gene" (using example file Additional file 2). 

Figure 5 shows an example plotting the magnitude 
of the effect size track in vertical format using the -b 
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Figure 6 Comparisons Across Groups. If the PheWAS analyses are stratified across multiple genetic ancestries or groups, it can be useful to 
view the similarities and differences of the significance of an association and direction of effect across groups. The output plot can be filtered by 
a single SNP (using -s SNP rs ID) and one or more groups (-r groupl, group2, . . .) by specifying an identifier for specific results in a column 
labeled "Groups" in the input file (example in Figure 6). Triangles point up for direction of effect that is positive, and point down for direction of 
effect that is negative. A different color represents each group. 



parameter, as well as plotting the sample size for all tests 
of association for each phenotype using the -A param- 
eter, with only points passing a p- value < 0.01 in blue. 

If the PheWAS analyses are stratified by population or 
genetic ancestry or other grouping, it can be useful to 
view the similarities or differences in the significance of 
an association and direction of effect across groups. The 
output plot can be filtered by a single SNP (using -s 
SNP rs ID) and one or more groups (-r groupl, group2, 
. . .) by specifying an identifier for specific results in a 
column labeled "Groups" in the input file. Figure 6 
shows results for the SNP rs673548, and African Ameri- 
cans (AA) and European Americans (EA), filtered by just 
"Allergy" phenotypes. Each group is represented by a dif- 
ferent color, and triangles point up for direction of gen- 
etic effect that is positive, and point down for direction 
of genetic effect that is negative. PheWAS -View recog- 
nizes the population abbreviations listed in Table 2 and 
used in Additional file 2. To use other populations or 

Table 2 PheWAS-View recognizes the population or 
genetic ancestry abbreviations listed below that can be 
used as "group" identifiers in the input file 



Label Color Description 



EA 


red 


European American 


AA 


blue 


African American 


H 


green 


Hawaiian 


API 


purple 


Asian Pacific Islander 


Al 


orange 


American Indian 



genetic ancestry or group abbreviations, an alternate 
group map file can be supplied by the user (-1 filename, 
txt) (example Additional file 4). 

Within PheWAS, significant phenotype-genotype 
results between SNPs may be due more to the relation- 
ship between phenotypes rather than the independent 
associations between a genetic variant and multiple phe- 
notypes (known as pleiotropy). For instance, if a series of 
cardiovascular disease-related measurements are highly 
correlated, the resultant SNP-phenotype associations 
may be very similar between all the cardiac disease phe- 
notypes due to the correlation between the phenotypes. 
PheWAS -View can help distinguish between correlated 
phenotypes and possible pleiotropy. If pairwise pheno- 
type correlations are exhaustively calculated and saved in 
tab-delimited matrix form, the absolute value of the cor- 
relation coefficients can be plotted using PheWAS -View 
by using the -C correlation file. Figure 7 A is an example 
PheWAS -View plot with the phenotypic correlation heat 
map, Figure 7B shows the same plot in vertical rather 
than horizontal format, where the cells of the correlation 
plot range from yellow to blue in the direction of 
decreasing absolute value of the correlations. Addi- 
tional file 5 is the example correlation matrix used for 
Figure 7A and B. 

Results can also be plotted with "expected" association 
results in blue, and novel associations plotted in purple 
by supplying a file of SNPs matched to individual pheno- 
types that are expected results (-x phenotype/SNP 
file), Figure 8, using Additional file 1 and Additional 
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(See figure on previous page.) 

Figure 7 Phenotypic Correlation Heatmap. A heatmap of correlations can be plotted with PheWAS-View if a matrix of all pairwise correlation 
coefficients between phenotypes is plotted. Figure 7 A is a PheWAS-View plot with the addition of a correlation heatmap. The cells of the 
correlation plot range from yellow to blue in the direction of decreasing absolute value of the correlations. Figure 7B shows the same plot in 
vertical rather than horizontal format after using the -a parameter. 
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Figure 8 Distinguishing "Expected" and "Novel" Associations Using Color. In this plot, novel associations with p < 0.01 are plotted in a 
different color (purple) from results that are more expected (blue). Using color to distinguish expected and more novel results facilitates 
prioritizing results to investigate further. 



PheWAS 



P hen B MO of intertrochenter region - gm/cm sq + 

1e-09 



Phen Cat - wheal length 



Phen Mononuclear percent + 
Phen German cockroach - wheal length + 
Phen Monocytes 
Phen Sum of deciduous or primary dft + 
Phen BMD ol trochanter region - gm/cm sq + 

Phen BMC of Wards triangle region - gm + 
Phen P amplitude, positive phase, lead V1 + 
Phen Bone minri density total regbn-gm/ctn sq + 




Phen Oral soft tissue: number of lesions + 

Phen German cockroach - flare width + 

Phen Natural Log Transformed - Serum insulin - 

Phen Bone mineral density femur neck -gm/cm sq + 
P hen Sum of permanent DMFS due to disease ■ 
P hen Bone area intertroc hanter region - cm sq + 

Phen Negative control - flare width + 
Phen Ragweed ■ flare length ■ 
Phen Radial puke + 



Phen Natural Log Transformed - Aspartate aminotransferase: SI + 



P hen Bone rn ineraf content of total region -gm + 



Phen BMC of trochanter region - gm + 



Figure 9 Sun Plot of Association Results for a Single SNP. PheWAS results plotted for a single SNP with p < 0.05 (-m p-value). The length 
of the line corresponds to the significance of the p-value, with the most significant result at the top ("noon") sweeping around clockwise. 
The p-value of the most significant result is listed to provide a sense of scale. The lines are red for p-values more significant than a threshold of 
p = 1x1 0" 3 (-p p-value). To add a "-" or "+" for direction of genetic effect for each phenotype use the flag -b. 
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file 6). Decisions about expected versus unexpected 
results are a study-by-study decision of the research- 
ers involved. One example would be to consider pre- 
viously reported SNP-phenotype associations as 
expected and previously unreported SNP-phenotype 
associations as novel Plotting the results with two 
different colors can be extremely useful to contrast 
results considered novel and those due to known 
genotype-phenotype relationships. 

PheWAS also provides an alternate way to visualize all 
significant results for a single SNP or phenotype through 
a "sun plot". Figure 9 shows a sun plot of all PheWAS 
results plotted for a single SNP with a p-value<0.05 
(using -m p-value to set the figure p-value cutoff). The 
length of the line corresponds to the significance of the 
p-value, where the more significant the p-value, the 
longer the line. The most significant result is at the top of 
the plot, with the p-value of the most significant result 
listed to provide a sense of scale. The remaining results 
sweep around clockwise. The lines are red for p-values 
more significant than a threshold of p = lxlO" 3 (using 
-p p-value to specify a threshold for lines being red or 
grey). To create this figure using Additional file 2, add- 
itional parameters were used: -S to create the sun plot, 
-s with SNPID to identify the specific SNP of interest for 
the sun plot, and -g to add the gene symbol to the plot. 
To add a "-" or for direction of effect for each associ- 
ation test, the -b flag was used with the Additional file 2, 
which contains direction of genetic effect information. If 
there is genetic ancestry or population information, using 
the -E flag will indicate the group for each of the asso- 
ciations. Other options for sun plots include plotting all 
results for a single gene at a p-value threshold (using -g 
gene name), or a single phenotype at a p-value threshold 
(using -P phenotype name) (not shown here). 

Conclusions 

The PheWAS approach provides a way to explore plei- 
otropy and the interrelationships between phenotypes, 
and as well as generate new hypotheses about the gen- 
etic architecture of complex traits. Visualizing complex 
PheWAS results with the various possible plots available 
within PheWAS -View provides a way to explore the data 
in a visual way, facilitating data analysis and interpret- 
ation. This software could be also be used for other 
phenotypically rich association studies such as expres- 
sion quantitative trait loci (eQTL) studies, studies that 
have high numbers of phenotypes due to the expression 
of multiple genes coupled with genotypic data. 

Availability and requirements 

Project name: PheWAS-View 

Project home page: http://ritchielab.psu.edu/ritchielab/ 
software 



Operating systems(s): Linux, Mac OS X, Windows 

Programming language: Ruby 

Other requirements: RMagick 

License: GNU General Public License 

Any restrictions to use by non-academics: The use 
of PheWAS -View is restricted to academic and non- 
profit users 

Additional files 



Additional file 1: The single input file used to produce the 
standard PheWAS-View plots of Figure 1, Figure 8, and Additional 
file 3. 

Additional file 2: The input file used to produce the PheWAS-View 
plots of Figures 3, 4, 5, 6, 7, and 9. Adding additional information 
columns allows for the use of additional features within PheWAS-View. 

Additional file 3: Figure SI Standard PheWAS-View Output 
Example Plot Sorted by P-value. A series of simulated PheWAS results 
plotted in PheWAS-View for a group of phenotypes. The y-axis presents 
-log 10 (p-value) of the tests of association for all SNPs for each 
phenotype, and the x-axis represents individual phenotypes. In this way 
all results for all SNPs are plotted for each phenotype. This plot was 
created using the same data as Figure 1, plotted sorting the input data 
file by p-value. 

Additional file 4: This is an example of the type of file that can be 
supplied by the user in order to use other population, genetic 
ancestry, or group abbreviations. 

Additional file 5: This is the example phenotype correlation matrix 
input file used for Figure 7A and Figure 7B. 

Additional file 6: The example file used to plot "expected" 
association results in blue, and novel associations in purple for 
Figure 8. 
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